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Abstract 

1-D scalar conservation laws with convex flux and Markov initial data 
are now known to yield a completely integrable Hamiltonian system. In 
this article, we rederive the analogue of Loitsiansky's invariant in hy- 
drodynamic turbulence from the perspective of integrable systems. Other 
relevant physical notions such as energy dissipation and spectrum are also 
discussed. 
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1 Introduction 

The inviscid Burgers equation 



d t u + d x \-jj =0, xeR,t>0, u(0,x)=u o (x) (1) 

with random initial data u has been extensively studied over the course of 
several decades. As a caricature of hydrodynamic turbulence, the earliest fun- 
damental results are due to Burgers himself, who studied (1) for uq a white noise 
process in space [7]. When considered with an added random forcing term, the 
equation becomes a model of random growth with deposition. Taken together, 
these form the basis of Burgers-KPZ turbulence. A broad overview of these 
areas can be found in [25] and the review article [1]. 

Many essential features of Burgers turbulence are preserved in the more 
general context of 1-D scalar conservation laws with strictly convex flux and 
random initial data. In recent work [21], we have shown the class of Markov 
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processes in space with only downward jumps is preserved by the deterministic 
dynamics. Although it is one of the simplest examples of a nontrivial flow on the 
space of probability measures, this model has surprisingly rich behavior. The 
evolution is given in terms of a Lax equation, and it has deep connections to 
kinetic theory, random matrices, and statistics. Moreover, it has been shown to 
be a completely integrable Hamiltonian system [19]. As an integrable system, 
its evolution is constrained by conserved quantities. In this article we examine 
how one particular conserved quantity of physical interest, well-known in the 
turbulence literature, arises in our setting. 



1.1 Entropy solutions to scalar conservation laws 

To begin, let us review some basic facts about the 1-D scalar conservation law 

d t u + d x f{u) = 0, iel,i>0 (2) 

with strictly convex flux / G C 1 and initial data u(x, 0) = u (x). Classical 
solutions found using the method of characteristics typically cease to exist after 
a finite time due to the formation of shocks. Instead, a suitable global-in-time 
solution to (2) is the entropy solution, derived via a variational principle for the 
action functional 

I(s;x,t) =U (s)+tf* (^r) ■ ( 3 ) 

Here, Uo(s) = Jq uo(r)dr is the initial potential and f*(s) = sup u6R {ws — f(u)} 
the Legendre transform of the flux. We assume that J7 has no upward jumps 
and that lim^i^^ I(s; x, t) = +oo always holds, ensuring that the infimum of 
/ is achieved. The inverse Lagrangian function a(x,t) is then defined by the 
Hopf-Lax formula: 

a(x, t) — arg + min I(s; x, t). (4) 

The + denotes that we choose a(x, t) to be the largest value where / achieves 
its minimum. In terms of infinitesimal particles, a gives the (rightmost) initial 
location of the particle at location x at time t. For fixed t > 0, a(x, t) is 
increasing in x. At points of continuity of a(-,t) the entropy solution is 

u(x,t) = (f')- 1 ( X ~ a t {X ' t) ). (5) 

Downward jumps in u correspond to shocks, and the velocity of a shock at x with 
left- and right-limits u± = u(x±,t) satisfies the Rankine-Hugoniot condition 

u( X ,t) = /("-)-/(*+) =; [f]u _ iU+ . ( 6 ) 

it_ — u + 
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1.2 Loitsiansky's invariant 

Now suppose that Uq(x) is a mean-zero velocity field which is statistically sta- 
tionary and ergodic in x. Then for each t > 0, the entropy solution u(x,t) 
remains stationary ergodic and ensemble averages are equivalent to spatial av- 
erages (which we denote by brackets). It has long been known that (1) admits a 
conserved quantity analogous to Loitsiansky's invariant in hydrodynamic tur- 
bulence [17]: 

f- oo 

Jo(t)= (u(x,t)u(x + h,t))dh. (7) 
Jo 

Jo is simply the integral of the two-point correlation function, and is essentially 
the value of the energy spectrum in the limit of low wavenumbers. In a series of 
articles [3, 4, 5, 6], Burgers provided three distinct derivations of the invariance 
of Jo- Upon doing so he investigated its implications on the evolution of a 
variety of shock statistics. Intriguingly, many of Burgers' calculations involve 
(unclosed) kinetic equations describing shock coalescence as in [11], and are 
quite prescient in that respect. 

It is easily seen that Jo is preserved for any 1-D scalar conservation law. A 
formal derivation is as follows. Denoting u(x) = u(x,t), one has that for any 
h > 0, 

d t (u(x)u(x + h)) = (u(x)d t u(x + h) + d t u(x)u(x + h)) 

= - (u(x)d x f(u(x + h)) + d x f(u(x))u(x + h)) . 

Therefore, the two-point correlation function 9(h) = (u(x)u(x+h)) itself satisfies 
a conservation law with flux r(h) = (u(x) (f(u(x + h)) — f(u(x — h)))}: 

d t e{h) + d hT (h) = o. (8) 

For Burgers' flux f(u) — u 2 /2, this is the analogue of the Karman-Howarth 
equation for three-dimensional isotropic turbulence in the vanishing viscosity 
limit. Assuming that the correlation length of the initial velocity field decays 
sufficiently fast, one has that r(h) — > as h — > oo and 

poo 

d t J (t) =dt 9(h)dh = 0. 
Jo 

If, in addition, u(x) is statistically invariant under the transformation x >-> —x, 
u i — y —u, the correlation flux takes the form 

r(h) = -2 (f(u(x))u(x + h)) . (9) 

When /(it) = u 2 /2, this implies that r = — -T/S3 with third-order structure func- 
tion 83(h) = ((u(x + h) — u(x)) 3 y This is particular to Burgers equation since 
t typically cannot be expressed in a simple manner using structure functions. 

Finally, we briefly remark on the physical relevance of the conserved quan- 
tity (7) for isotropic Navier-Stokes turbulence in the infinite domain. In this 
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context, it is widely believed that certain universal scaling laws, such as a law 
of energy decay should hold for the statistics of the velocity field. Dimensional 
analysis dictates that the mean energy dissipation should be proportional to 
U 3 /L, where U(t) is the root-mean square velocity and L(t) is the integral 
length scale. A self-similarity ansatz for the velocity correlation function [24] 
relates the scaling of U(t) and L(t) — however, an additional relation is needed 
to obtain the scaling exponents themselves. The invariance of Jo, commonly 
referred to as the "permanence of large eddies," was utilized by Kolmogorov 
[16] to obtain a power law for the energy decay. While this was later shown 
to hold only under particular assumptions on the initial spectrum, it remains a 
fundamental result in the theory. In the simpler setting of Burgers turbulence, 
many of the same issues have been thoroughly investigated and settled. We 
direct the interested reader to [1, 14, 15] for a comprehensive discussion. 

1.3 Lax equation and integrability of shock clustering 

To start, we take uq{x) to be a stationary Markov process with only downward 
jumps. In [21], it was proven that the entropy solution to (2) then satisfies a 
closure property: for each fixed t > 0, u{x,i) is a stationary Markov process in 
x with only downward jumps. This was also shown to hold if the initial data 
is the derivative of a stationary Levy process in x with only downward jumps 
(e.g., white noise). A formal derivation of this closure for Burgers equation was 
obtained in [10]. 

Denote by p(dy,t) and {qh{y,dz,t)} h>0 the stationary and transition mea- 
sures in x of the process u(x, t). The evolution of the statistics can be stated in 
terms of the generator of the solution process. First, define the 1- and 2-point 
operators through their action on appropriate test functions ip: 

V{t)ip{y)= [ p(dy,t)ip(y), Qh(t)tp(y) = [ q h {y,dz,t)ip{z). (10) 
Jr Jr 

Consider the generator A(t)ip = lim^o ji (Qh(t)f — <p) of the Markov semigroup 
in x (assuming the Feller property holds). For each fixed t > 0, it has the form 

A(t)<p(y) = b{y, t)v'(y) + f n{y, dz, t)(<p(z) - <p(y)), (11) 

Jr 

where b(y,t) is the drift coefficient of the process and n(y,dz,t) is the jump 
measure. Now we define an operator B(t) which, for each fixed x, serves as a 
'generator' of the solution in t: 

B(t)<p(y) = -b(y,t)f'{v)<p'{y) - f n(y,dz,t)[f] y , z (<p(z) - <p(u)). (12) 

Jr 

B(t) truly is the generator of a Markov semigroup if / is decreasing. As shown 
in [21], the evolution of the 1- and 2-point operators then takes the form 

d t V = VB, d t Q h = [Q h ,B] (13) 
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where square brackets denote the commutator. In terms of generators, this is 
the Lax equation 

d t A=[A,B]. (14) 

We have shown that (14) is equivalent to a kinetic equation for the shock statis- 
tics. These equations (remarkably!) admit explicit solutions when uq is a Brow- 
nian motion [2, 8, 9, 22] or a white noise [7, 12, 13]. In addition, there are many 
intriguing links to random matrices and integrable systems. This is discussed 
at length in [21]. 

2 The Manakov relation and an invariant of the 
Lax equation 

While the Lax equation yields some conserved quantities, these are not sufficient 
to fully describe the evolution of the system. One can, in addition, verify that 
the operators A and B satisfy the Manakov relation 

[A, AT] - [M,B] = 0. (15) 

with multiplication operators 

M<p(y) = y<p(y), Nip{y) = f(y)<p{y). (16) 

Equation (15) is rigid in that it only holds with multiplication operators M. and 
M: if we let M^,<p(y) = ip(y)ip(y) for any ip, then [A, M^ A ] - [M^ B ,B] = if 
and only if iPa{v) = f(y) and ip B (y) = V- 

The Manakov relation yields the necessary additional integrals by allowing 
a spectral parameter to be introduced as in [18]. That is, by (14) and (15), 

d t (A-nM) = [A- iiM,B + nN], fieC. (17) 

If one considers discretizations of A and B by TV x TV matrices, this implies that 
the spectral curve 

r = { (A, n) e C 2 1 det(A - Aid - pM) = 0} . (18) 

remains unchanged in time. The conserved quantities are then given by the co- 
efficients of the characteristic polynomial. It has recently been shown by Menon 
[19] that this discretized system, the Markov TV-wave model, is a completely 
integrable Hamiltonian system. While the argument above shows that the spec- 
trum is invariant in the discrete setting, this does not remain true as TV — > oo 
since the continuous spectrum may evolve. Therefore, it is an interesting and 
open problem to find conserved quantities that survive in the continuum limit, 
and to determine which of these are finite. Since the discrete model does not 
allow for nontrivial statistically stationary solutions unless one takes TV — > oo, 
we note that the invariant Jo can only appear in this limit. 
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Let us now state our main result. In what follows we will assume that p{dy, t) 
and qh{y, dz, t) decay fast enough for the appropriate integrals to converge. For 
stationary Markov solutions to (2), we find that the invariance of 

Jo= f E[u(x,t)u(x + h,t)]dh = f ( f p{dy,t)y f q h {y, dz,t)z) dh (19) 
Jo Jo \Jr Jr J 

is equivalent to the identity 

VA(.f(y)Q h y) = 0. (20) 

Since u is stationary in x, this is the weak form of the forward equation = A^p 
with test function f(y)Qh,y- More generally, fix n € N, denote hi — Xi — a^-i 
and Qi = Qh t for i = 1, . . . , n — 1, and let h = (hi,..., h n -i)- We demon- 
strate that the Manakov relation (15) implies the n-point function (9(h, t) = 
VyQ\y ■ ■ ■ Qn-iV satisfies the conservation law 

d t 0(h) + V • T(h) = (21) 

with flux Ti(h, t) = VyQxy ■ ■ ■ Qi-i (/(»)&!/ - yQif(y))Qi+iV ■ • ■ Qn-iV- lute- 
gration with respect to h shows that 

Jo n) = / E[u(x ,t)u(xi,t)---u(x n -i,t)]dxi---dx n -i (22) 

Jjo<Il<"<I„-l 

is conserved by the evolution. 

We note that these invariants can also be arrived at by the formal argument 
in the introduction, as the essential properties are the stationarity of the velocity 
field and sufficient decay in the correlation length scale. The purpose of this 
article, however, is not to prove a new result but to provide a novel perspective 
on the matter from the viewpoint of integrable systems. 



2.1 Derivation of main result 

For simplicity we will suppress t in the notation whenever possible. We begin by 
showing the equivalence of (19) and (20). To do this, we need only demonstrate 
that the conservation law for the one-point function 

d t (TyQ h y) = -d h (V {-f(y)Q h y + yQ h f(y)}) (23) 

holds. Recall that dhQh = AQh = QhA. Then (23) is equivalent to 

V {.f{y)AQ h y - yQ h Af(y) - ByQ h y + y[Q h , B]y} = 0. (24) 

A short computation gives that 

f(y)AQ h y - yQ h Af(y) =b(y)f(y)(d y Q h y) - yQ h (b(y)f(y)) 

+ f(y) [ n(y,dz)(Q h z-Q h y) 
Jm 

-vQh(f n(y,dz)(f{z)-f(y))) . 
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Similarly, 

ByQ h v + y[Q h , B]y = - b(y)f(y)Q h y - yQ h (b(y)f'(y)) 
- [ n(y,dz)(f(z)-f(y))Q h z 

-yQh Qf nfod*)(/(*)-/(l/))) ■ 

Substituting these into (24) and collecting terms yields (20). 

Similar calculations hold for the n-point function and become clearer upon 
explicit use of the Manakov relation (15). Using that constants are in the 
nullspace of A and £>, equations (14) and (16) imply that 

dt{VyQxy-Q n -iy) 

= VByQty ■ ■ ■ Q n ^y + ^ PyQ iy ■ ■ ■ Qi-iy[Qi, B]yQ l+1 y ■ ■ • Q„_ l2/ 
i=i 

n-l 

= V[B, M]Qiy ■ ■ ■ Q n -xy + £ VyQ iy ■ ■ ■ Q l [B, M] Q l+1 y ■ ■ ■ Q n - lV . 

i=l 

Substituting (15), the previous expression becomes 

n-l 

V[M, A]Qiy ■ ■ ■ Qn-xy + ^ VyQ iy ■ ■ ■ Qi-iyQi[Af, A]Q t+1 y ■ ■ ■ Q n - X y 

i=l 

n-l 

= ^ J d hi {VyQvy ■ Q l -i{f{y)Q l y - yQif(y))Q i+ iy ■ ■ ■ Q n -iy}, 
i=i 

which is equivalent to (21). 

We remark that the above argument clarifies why quantities not of the form 
(22), such as jE[u 2 (x)u(x + h)]dh, are typically not conserved in time: 

/>oc />oo 

dt Vy 2 Q 2 ydh 2 = d t VyQ iy Q 2 y\ hi=0 dh 2 
Jo Jo 

t>OC 

= / V{f{y)Ay - yAf(y))Q 2 ydh 2 ? 
Jo 

except when y and f(y) are in the nullspace of A. 

3 Exact solutions, energy, and spectra 
3.1 White-noise initial data 

The invariant Jo is expressed in terms of the 1- and 2-point functions for solu- 
tions which retain the Markov property in space. In certain special cases, its 
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value can be explicitly obtained. For example, it is straightforward to compute 
Jo when uo(x) is a stationary Ornstein-Uhlcnbeck (mean- reverting) process: 

duo(x) = —(3u (x)dx + jdB X} /3, 7 > 0, (25) 

where B x is Brownian motion. Since u has zero mean, 

2 

E[u {x)u (x + h)] = Cov(u (x),u (x + h)) = ^e~ 0h . (26) 

Therefore, J = 7 2 /(2/3 2 ). 

By taking appropriate limits we obtain Jo in the case of a white noise initial 
velocity field for any flux /. That is, suppose the initial potential is Uq(x) = uB x 
with a > the strength of the noise and B x a two-sided Brownian motion pinned 
at the origin. This corresponds to letting 7 — afi and taking /3 — > 00 in (25), so 
that Jo = cr 2 /2. The solution to Burgers equation with white noise initial data 
was obtained explicitly in terms of a Painleve transcendent by Groeneboom [13], 
and later rediscovered by Frachebourg and Martin [12]. The above computation 
thus determines the value of an integral of the form (19) with stationary and 
transition measures as in [12, 13] . For a more thorough examination of physically 
relevant quantities we refer the reader to [23]. 



3.2 Energy dissipation and power spectrum 

Jo is closely related to other quantities of physical interest, including the mean 
energy dissipation per unit interval and power spectrum. We now discuss some 
basic facts regarding these quantities for stationary Markov solutions to (2). 

The mean dissipation in an interval I is computed as follows. By Fubini's 
theorem and equation (13), 



d 4 E 



w\IM x,t)i 



dx 



\dtVy 2 - \VB V 2 . 



Let h(y) = wf'{w)dw be the entropy flux associated to y 2 /2. Then 



l -By 2 = -b(y)yf'(y) - f n(y,dz){f] y . 



= ~Ah{y) 



_ - VI 

2 2 

2 2 

Applying V and using that VAh — 0, we obtain 
1 



n(y,dz) \ [f] Vt 



(Hz) - h{y)) 



d t E 



-u(x,t) dx 



= - I P(dy) I n(y,dz)\hz + y)(f(z)-f(y)) 

JR JR L z 



(27) 



(28) 



wf'(w)dw \ . (29) 
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This can also be arrived at by considering traveling wave solutions to (2) as in 
[20]. In particular, the energy dissipated at any instant by a shock connecting 
the states u_ > u + is 



1 r u - 

-(«+ + «_)(/(«+)-/(«_)) + J wf> 



(w)dw 



= \{u- - «+)(/(«_) - /(«+)) - jf /Hdw > 0. 

Summing up over the expected number of shocks from w_ to u + per unit interval 
implies (29). 

Next, consider the power spectral density of the velocity field u(x,t): 

E(k,t) = -!- / e- lkh K[u(x,t)u(x + h,t)]dh, keR,t>0. (30) 

Again, we drop t from the notation. For stationary Markov solutions, the power 
spectrum is well-defined and can be given in terms of the A-resolvent of the 
transition semigroup (Qh)h>o- With ip a test function, let 

/>oo 

KMy) = / e- Xh Q h ip(y)dh, Re A > 0. (31) 
Jo 

Extending 1Z\ to A = ik with ^ k G R, we have that 

E(k, t) = -Re {V{t)yTl lk {t)y) . (32) 

7T 

The Laplace transform of (23) with respect to h is therefore seen to be equivalent 
to an evolution equation for the power spectrum: 

dtE = --Re (P{f(y)1l ik y - y1l ik f(y)}) . (33) 

To obtain IZik in (32), one first needs to solve the Lax equation (14). Finally, 
we note that Jo = lim^o irE(k, t) is simply a constant multiple of the spectral 
density in the low-wavenumbcr limit. 

Let us contrast the situation considered here with that of Burgers equation 
with one-sided Levy initial data as in [20]. In the latter case, it was shown 
that there exist solutions with finite energy and infinite dissipation per unit 
interval due to an influx of energy through the boundaries. Here, the stationarity 
of the solution forbids such energy fluxes. It was also demonstrated in [20] 
that the traditional notion of power spectrum is irrelevant when considering 
solutions with independent increments, as processes with the same energy have 
indistinguishable spectra. A useful notion of spectrum was instead given by the 
Fourier-Laplace transform of process paths without averaging. In contrast, here 
we have that the power spectrum is nontrivial and can be computed via the 
resolvent of the solution process. 
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